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Abstract 

The Gamow- Teller transition operator is written as a polynomial in the dipole proton-neutron 
and quadrupole charge conserving QRPA boson operators, using the prescription of the boson 
expansion technique of Belyaev-Zelevinski type. Then, the 2v(3(3 process ending on the first 2 + 
state in the daughter nucleus is allowed via one, two and three boson states describing the odd-odd 
intermediate nucleus. The approach uses a single particle basis which is obtained by projecting out 
the good angular momentum from an orthogonal set of deformed functions. The basis for mother 
and daughter nuclei have different deformations. The GT transition amplitude as well as the half 
lives were calculated for ten transitions. Results are compared with the available data as well as 
with some predictions obtained with other methods. 

PACS numbers: 23.40.Hc,21.10.Tg,21.60.Jz,13.10.+q,27.50.+e 
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One of the most exciting subject of nuclear physics is that of double beta decay. The 
interest is generated by the fact that in order to describe quantitatively the decay rate one 
has to treat consistently the neutrino properties as well as the nuclear structure features. 
The process may take place in two distinct ways: a) by a 2vf3f3 where decay the initial nuclear 
system, the mother nucleus, is transformed in the final stable nuclear system, usually called 
the daughter nucleus, two electrons and two anti-neutrinos; b) by the 0i//3/3 process where 
the final state does not involve any neutrino. The latter decay mode is especially interesting 
since one hopes that its discovery might provide a definite answer to the question whether 
the neutrino is a Majorana or a Dirac particle. The contributions over several decades have 
been reviewed by many authors, [l, 2, 3), 0, Q]. 

Although none of the double beta emitters is a spherical nucleus most formalisms use a 
single particle spherical basis. 

In the middle of 90 's we treated the 2vf}(3 process in a pnQRPA formalism using a 
projected spherical single particle basis which resulted in having a unified description of the 
process for spherical and deformed nuclei [7, 8]. Recently the single particle basis [j], [h| has 
)een improved by accounting for the volume conservation while the mean field is deformed 
111 . \\\ . The improved basis has been used for describing quantitatively the double beta 
decay rates as well as the corresponding half lives Q] • The results were compared with 
the available data as well as with the predictions of other formalisms. The manners in which 
the physical observable is influenced by the nuclear deformations of mother and daughter 
nuclei are in detail commented. Two features of the deformed basis are essential: a) the 
single particle energy levels do not exhibit any gap; b) the pairing properties of the deformed 
system are different from those of spherical system. These two aspects of the deformed nuclei 
affect the overlap matrix of the pnQRPA states of mother and daughter nuclei. Moreover, 
considering the Gamow- Teller (GT) transition operator in the single particle-space generated 
by the deformed mean-field, one obtains an inherent renormalization with respect to the one 
acting in a spherical basis. 

In Ref. ( 15() we studied the higher pnQRPA effects on the GT transition amplitude, by 
means of the boson expansion technique for a spherical single particle basis. Considering 
higher order boson expansion terms in the transition operator, significant corrections to 
the GT transition amplitude are obtained especially when the strength of the two body 
particle-particle (pp) interaction approaches its critical value where the lowest dipole energy 
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is vanishing. As we showed in the quoted reference, there are transitions which are forbidden 
at the pnQRPA level but allowed once the higher pnQRPA corrections are included. An 
example of this type is the 2v(3f3 decay leaving the daughter nucleus in a collective excited 
state 2 + . The electrons resulting in this process can be distinguished from the ones associated 
to the ground to ground transition by measu ring , in coincidence, the gamma rays due to the 



transition 2 + — > + in the daughter nucleus 16]. 



The aim of this letter is to study the double beta decay + — > 2 + where + is the ground 
state of the emitter while 2 + is a single quadrupole phonon state describing the daughter 
nucleus. The procedure is the boson expansion method as formulated in our previous paper 
[15I but using a projected spherical single particle basis. 

In order to fix the necessary notations and to be self-contained, in the present work we 
describe briefly the main ideas underlying the construction of the projected single particle 
basis. 

The single particle mean field is determined by a particle-core Hamiltonian: 

H = H m + H eore -Mdy £ E <W» (0.1) 

A=0,2 ~\<fj,<\ 

where H sm denotes the spherical shell model Hamiltonian while H core is a harmonic 
quadrupole boson (6+) Hamiltonian associated to a phenomenological core. The interac- 
tion of the two subsystems is accounted for by the third term of the above equation, written 
in terms of the shape coordinates «oo, otin- The quadrupole shape coordinates are related to 
the quadrupole boson operators by the canonical transformation: 

a2M= fc v 7 2 ( ^ + HV/,) ' (0,2) 

where k is an arbitrary C number. The monopole shape coordinate is determined from the 
volume conservation condition. In the quantized form, the result is: 

1 



ctoo 



(0.3) 



Averaging H on the eigenstates of H sm , hereafter denoted by \nljm), one obtains a deformed 
boson Hamiltonian whose ground state is, in the harmonic limit, described by a coherent 
state 

* g = exp[d(b% -b 2O )]\0) b , (0.4) 
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with |0)& standing for the vacuum state of the boson operators and d a real parameter which 
simulates the nuclear deformation. On the other hand, the average of H on \1/ g is similar 
to the Nilsson Hamiltonian [l^]. Due to these properties, it is expected that the best trial 
functions to generate a spherical basis are: 

K% = \rdjm)9 a . (0.5) 

The projected states are obtained by acting on these deformed states with the projection 
operator 

Pmk = J Di IK *(Q)R(Q)dQ. (0.6) 

The subset of projected states : 

KfM = KjPLiimm,] = Ki^TM' (°- 7 ) 

are orthogonal with the normalization factor denoted by A/Jjy. 

Although the projected states are associated to the particle-core system, they can be 
used as a single particle basis. Indeed, when a matrix element of a particle like operator is 
calculated, the integration on the core collective coordinates is performed first, which results 
in obtaining a final factorized expression: one factor carries the dependence on deformation 
and one is a spherical shell model matrix element. 

The single particle energies are approximated by the average of the particle-core Hamil- 
tonian H' = H — H core on the projected spherical states defined by Eq. (I0.7|) : 

^HOOI^Wu (<*)>■ (o-8) 

The off-diagonal matrix elements of H' is ignored at this level. Their contribution is however 
considered when the residual interaction is studied. 

As shown in Ref. |9[], the dependence of the new single particle energies on deformation is 



similar to that shown by the Nilsson model [17[ • The quantum numbers in the two schemes 
are however different. Indeed, here we generate from each j a multiplet of (2j + 1) states 
distinguished by the quantum number I, which plays the role of the Nilsson quantum number 
Q and runs from 1/2 to j and moreover the energies corresponding to the quantum numbers 
K and -K are equal to each other. On the other hand, for a given I there are 2/+1 degenerate 
sub-states while the Nilsson states are only double degenerate. As explained in Ref.[9|, the 
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redundancy problem can be solved by changing the normalization of the model functions: 

($™|$™> = i £($i A/ |$i M > = 2 - (o-9) 

M 

Due to this weighting factor the particle density function is providing the consistency result 
that the number of particles which can be distributed on the (21+1) sub-states is at most 
2, which agrees with the Nilsson model. Here a stands for the set of shell model quantum 
numbers nlj. Due to this normalization, the states & a M used to calculate the matrix elements 



of a given operator should be multiplied with the weighting factor J 2/ (21 + 1). 



Finally, we recall a fundamental result, obtained in Ref . [12| . concerning the product of 
two projected states, which comprises a product of two core components. Therein we have 
proved that the matrix elements of a two body interaction corresponding to the present 
scheme are very close to the matrix elements corresponding to spherical states projected from 
a deformed state consisting of two spherical single particle states times a single collective 
core wave function. The small discrepancies of the two types of matrix elements could be 
washed out by using slightly different strengths for the two body interaction in the two 
methods. 

As we already stated, in the present work we are interested to describe the Gamow- Teller 
two neutrino double beta decay of an even-even deformed nucleus. In our treatment the 
Fermi transitions, contributing about 20% to the total rate, and the "forbidden" transitions 
are ignored, which is a reasonable approximation for the two neutrino double beta decay in 
medium and heavy nuclei. The 2vf3(3 process is conceived as two successive single (3~ tran- 
sitions. The first transition connects the ground state of the mother nucleus to a magnetic 
dipole state 1 + of the intermediate odd-odd nucleus which subsequently decays to the first 
state 2 + of the daughter nucleus. The second leg of the transition is forbidden within the 
pnQRPA approach but non- vanishing within a higher pnQRPA approach [ijj]. The states, 
involved in the 2vf3f3 process are described by the following many body Hamiltonian: 

2 G 

H = 2J+~1 ^ TaI ~ ^ Ta ) CraIMCTaIM ~ X/ ~^^toliPtcxV 

- £ X ry Q T Ql. (0.10) 

T,r'=p,n 

The operator c' TOtIM (c TOl iM) creates (annihilates) a particle of type r (=p,n) in the state 
$^ M , when acting on the vacuum state |0). In order to simplify the notations, hereafter 



the set of quantum numbers a(— nlj) will be omitted. The two body interaction consists 
of three terms, the pairing, the dipole-dipole particle hole (ph) and the particle-particle 
(pp) interactions. The corresponding strengths are denoted by G T , x, Xi, respectively. All 
of them are separable interactions, with the factors defined by the following expressions: 
P t / = E 2 j j 



21 + \ t1mC tLM 



M 



M,M' * * 
MM' 1 1 



QS? = E4 T) (4c fe ) <&> = J-^—(h\\r 2 Y 2 \\I k ). (0.11) 

i,k ' V ~r 

The remaining operators from Eq. fl0.10l) can be obtained from the above defined operators 
by hermitian conjugation. 

The one body term and the pairing interaction terms are treated first through the 
standard BCS formalism and consequently replaced by the quasiparticle one body term 
YLriM E T Q>l-iM a TiM- in terms of quasiparticle creation (a rIM ) and annihilation (a r /jvf) oper- 
ators, related to the particle operators by means of the Bogoliubov-Valatin transformation, 
the two body interaction terms, involved in the model Hamiltonian, can be expressed just 
by replacing the operators (3.2) by their quasiparticle images. Thus, the Hamiltonian terms 
describing the quasiparticle correlations become a quadratic expression in the dipole and 
quadrupole two quasiparticles and quasiparticle density operators: 

m p ,m n 

Bl^ipn) = E ^m p -m„ ^ a P I pmp a nI n rn n (~) " "i 
rn p ,m n 

"^2/n( rr ) = £ ^mjm^, ^ a il T m T a i'I T ,m, T ,y 
m T ,m T i 

BUtt') = £ Ct^Ai^r'i^A-Y^'. (0.12) 

m T ,m T i 

Since the pnQRPA treatment of the dipole-dipole interaction in the particle-hole (ph) 
and pp channels run in an identical way as in our previous publications [l3, [ijj], here we 
don't give any detail about building the dipole proton-neutron phonon operator : 

rl = £[*i(*MJm(*) - (o.i3) 



We just mention that the amplitude are determined by the pnQRPA equations and the 
normalization condition. 

The charge conserving QRPA bosons 

= £[X 2 (£04#) - li(A0A,-#)(-n k = (jp,p'), (n,n'), (0.14) 
k 

are determined by the QRPA equations associated to the matrices: 

A TT ^t, i'k') = 5 TT ,5 tl >5 kk ,(Ei + ED - x TT , (4 r) 4 T) ) (g£eS) , 

B^Cifc; = -X TT , (qft&l) ,<<*,<'< k'. (0.15) 



+ (0,6) 

'1 + kk 



where 



In order to distinguish between the phonon operators acting in the RPA space associated to 
the mother and daughter nuclei respectively, one needs an additional index. Also, an index 
labeling the solutions of the RPA equations is necessary. Thus, the two kinds of bosons will 
be denoted by: 

rU?*:), j = ij; k= 1,2,. ..N^;Tl(jk), j = i,f- fc = l,2,..JV< 2 >. (0.17) 

Acting with r{ M (zfc) and T\ Jfk) on the vacuum states |0)j and |0)/ respectively, one obtains 
two sets of non-orthogonal states describing the intermediate odd-odd nucleus. By contrast, 
the states r^i/c) |0)j and r^C/^) |0) / describe different nuclei, namely the initial and final 
ones participating in the process of 2v(3(3 decay. The mentioned indices are however omitted 
whenever their presence is not necessary. 

Within the boson expansion formalism the transition GT operators are written as 
polynomial expansion in terms of the QRPA boson operators with the expansion coeffi- 
cients determined such that the mutual commutation relations of the constituent operators 
A\jjm), Ai tl (pri), Bl fl {pn), Bi tl {pn) be preserved in each order of approximation 18]. One 
arrives at the expressions: 

Al(hJn) = £ K;%Wn)rl#i) +A k ; 1 \j P jn)T 1 ^(k 1 )(-r»} 

ki 

+ £ UtZiiwn) [(rUks^M) r\(ki)] + A$f£ kl (j p j n ) [(r^r^))^)] 

k lt k2,k 3 ;l=0,2 1 ' » 
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+ E UZIUhJn) [r!(*o (rt(fc 2 )r 2 (fc 3 )) J 1 + 4gl(^) [(rj(fe)r 2 (fe)) r^) 

fcl,fc 2 ,*3;J=0,2 1 M ' 
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= E feCvn) frKfcOrJc^)] +Bgg(j l j„) [r^ry^)] 



fcifc 2 



+ sggfCwO [r{(*i)r 2 (fe) 



4, 



where the expansion coefficients are those given in Ref. [15j . If the energy carried by leptons 
in the intermediate state is approximated by the sum of the rest energy of the emitted 
electron and half the Q- value of the double beta decay process 

1 



AE = m e c 2 + X -Q^ 2) , 



the reciprocal value of the 2vf3f3 half life can be factorized as: 



(0.19) 



2+r 



F 2 |M^| 2 , 



(0.20) 



where _F 2 is the Fermi integral which characterizes the phase space of the process while the 
second factor is the GT transition amplitude which, in the second order of perturbation 
theory, has the expression: 

i(0H/3 + ||fc, m )j l (k,m\k',m') f f (k',m'\\P + \\2f) f 

k,m 



M, 



(02) 



GT 



v^E 



(Ek,m + AE 2 ) 



(0.21) 



Here AE 2 = AE + Ei+ , with E%+ standing for the experimental energy for the first state 1 + . 
The intermediate states \k,m) are k-boson states with k — 1,2,3 labeled by the index m, 
specifying the spin and the ordering label of the RPA roots. Inserting the boson expansions 
from Eq. fl0.18l) into the expression of the f3 + transition operator one can check that the 
following non- vanishing factors, at numerator, show up: 

^olliMi^OHi, i fel Mi, i fc2 ||ri(/, fc 2 )r 2 (/, i)|[i,2 1 > / , 

j <0||r 1 (^A;i)r 2 (i,A; 2 )||2,l fcl 2 fca ) i/ <2,l il 2 1 ||rl(/,j 1 )||l,2 1 ) / , 

2fc 2 2fc 3 )i/{3, l 3 - 1 2j 2 2 1 ||r i (/,ji)r 2 (/, j 2 )||l, 2 1 ) f , 
J (0||r 1 (z,A; 1 )r 2 (2,A; 2 )||2,l fel 2 fc2 ) j/ (2,l, 1 2 j2 ||r t 1 (/,j 1 )r 2 (/,j 2 )r 2 (/, 1)111,2!)/. (0.22) 



The term E^ m from the denominator of Eq. (10.211) is the average of the energies of the 
mother and daughter states \k,m) normalized to the average energy of the first pnQRPA 
states 1 + in the initial and final nuclei. Calculations were performed for the following 10 



(0.18) 
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double beta emitters: 48 Ca, 96 Zr, 100 Mo, 104 Ru, 110 Pd, 116 Cd, 128 Te, 130 Te, 134 Xe, 136 Xe. Since 
the single particle space, the pairing interaction treatment, and the pnQRPA description 
of the dipole states describing the intermediate odd-odd nuclei used in the present paper 



are identical with those from Refs. 



14| for ground to ground transition, we don't present 



them again. The strength of the QQ interaction was fixed by requiring that the first root of 
the QRPA equation for the quadrupole charge conserving boson is close to the experimental 
energy of the first 2 + state. The results of the fitting procedure are given in Table I. 

Having the RPA states defined, the GT amplitude has been calculated by means of 
Eq. ( 10.211) . while the half life with Eq. ( 10.201) . The Fermi integral for the transition + — > 2 + 
was computed by using the analytical result given in Ref.(|4|). 

The final results are collected in Table II. Therein one may find also the available experi- 
mental data as well as some theoretical results obtained with other approaches. One notices 
that the half life is influenced by both the phase space integral (through the Q-value) and 
the single particle properties which determine the transition amplitude. Indeed, for 128 Te 
and 134 Xe the small Q-value causes a very large half life, while in 48 Ca the opposite situation 
is met. By contrary the Q value of 110 Pd is about the same as for 76 Ge but, due to the 
specific single particle and pairing properties of the orbits participating coherently to the 
process, the half life for the former case is more than three orders of magnitude less than 



in the later situation. Since in Ref. 



19| a spherical single particle basis and a renormalized 



pnQRPA approach are used, the transition matrix elements are larger and the half lives 
shorter than in our calculations, although the higher RPA formalism in the quoted reference 

n n 

is similar [21] to the one proposed by us in Ref. [15] and used in the present work. It is 
worth mentioning the good agreement between our prediction for 100 Mo and that of Ref. [^(J 
obtained with a deformed SU(3) single particle basis. 

It is worth mentioning that the double beta transitions to excited states have been consid- 
ered by several authors in the past, but the calculations emphasized the role of the transition 
operator and some specific selection rules. Many of calculations regarded the neutrinoless 
process. Thus, in Ref. 22J it was shown that the neutrinoless transition to the excited + for 
medium heavy nuclei might be characterized by matrix element which are larger than that of 
ground to ground transition and that happens since in the first transition, the change of the 
K quantum number is less. In Ref. [23!] it has been stated that the + — > 2 + matrix element 
depends on the left-right current coupling and not on the neutrino mass. However according 
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Nucleus 


E e 9 x + p - [keV] 


E& [keV] 


b A X w [keV] 


48 Ca 


983 


983 


0.0713 


48 Ti 


983.52 


979.02 


42.8 


76 Ge 


562.93 


558.88 


50.8 


76 Se 


559.10 


558.87 


65.2 


96 Zr 


1750.49 


1465.6 


2. 


96 Mo 


778.24 


776.8 


38.1 


100 Mo 


535.57 


534.4 


31.5 


100 Ru 


539.5 


536.1 


19.7 


104 Ru 


358.03 


358.45 


29.8 


104p d 


555.81 


561.83 


20.9 


no pd 


373.8 


370.45 


44.65 


110 Cd 


657.76 


662.8 


25.1 


116 Cd 


513.49 


514.5 


30.5 


116 Sn 


1293.56 


1179.16 


7.0 


128 Te 


743.22 


746.12 


12.12 


128 Xe 


442.91 


449.58 


19.43 


iao Te 


839.49 


831.03 


12.12 


130 Xe 


536.07 


534.2 


17.28 


134 Xe 


847.04 


841.75 


20.0 


134 Ba 


604.72 


607.98 


17.56 


136 Xe 


1313.027 


1314.9 


16.37 


136 Ba 


818.49 


810.3 


14.82 



TABLE I: The experimental and calculated energies for the first 2 + states in mother and daughter 
nuclei are given. The strength parameter of the quadrupole-quadrupole interaction was fixed so that 



the experimental energies are reproduced. In our calculations we considered X pp = X 



X, 



pu ■ 



The oscillator length is denoted by b = (Ti/Moj) 1 / 2 . 

to the calculations of Haxton et a/[2l], the matrix element is suppressed and therefore it is 
not possible to extract a meaningful parameter for the left-right coupling. Although the 
transition operator might have a complex structure, many calculations have been performed 
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Nucleus 


o 2+ 

RR 
^ PP 


AE 2 [MeV] 


\M { rr 2) \[MeV- 1 ] 


T^[yr] 




units of m e c 2 






present 


Exp. 


Suhonen [19] 


48 Ca 


6.432 


2.473 


0.901-10- 3 


1.72-10 24 






76 Ge 


2.894 


1.295 


0.558-10~ 3 


5.75-10 28 


>1.M0 21 


1.0-10 26 


96 Zr 


5.033 


2.913 


0.834-10" 3 


2.27-10 25 


> 7.9-10 19 


4.8-10 21 


100 Mo 


4.874 


1.756 


0.136-10- 2 


1.21-10 25 


>1.6-10 21 


3.9-10 24 
a )2.5-10 25 
6 )l.2-10 26 


104 Ru 


1.456 


0.883 


0.028 


6.2-10 28 






no pd 


2.646 


1.182 


0.050 


1.48-10 25 






116 Cd 


2.967 


1.269 


0. 507-10~ 2 


3.4-10 26 


>2.3-10 21 


1.1-10 24 


128 Te 


0.836 


1.305 


0.229-10" 2 


4.7-10 33 


>4.7-10 21 


1.6-10 30 


i30 Te 


3.902 


2.358 


0.620-10~ 3 


6.94-10 26 


^4.5-10 21 


2.7-10 23 


134 Xe 


0.460 


0.806 


0.621-10~ 2 


5.29-10 35 






136 Xe 


3.251 


1.518 


0.249-10" 2 


3.88-10 26 




2.0-10 24 



TABLE II: The results for the GT transition amplitudes as well as for the half lives of the double 
beta decay + — > 2 + are given. Also the Q values of the transitions are given in units of m e c 2 . AE 2 
is the energy shift defined in the text. For comparison, we give also the available experimental 
results as well as some theoretical predictions obtained with other formalisms. For 100 Mo we 
mention the result of Ref. ( 



20]) obtained with an SU(3) deformed single particle basis a ) and with 



a spherical basis b \ 



with the approximate interaction [cr(l) x cr(2)] A=2 t + (l)i + (2) in order to test some selection 
rules. Thus, this interaction was used in Ref. [241] for the transition + — > 2 + of 48 Ca, us- 
ing a single j calculation. It has been proved that the matrix element for this transition 
is suppressed due to the signature selection rules. The transition to 0^ was examined for 
A = 76, 82, 100, 136 nuclei by assuming light and heavy Majorana neutrino exchange mech- 



anism and triliniar R-parity contribution. Hig 



the nuclear matrix elements were included 



ler RPA as well as renormalization effects for 



25]. 



Here we show that the transition + — ► 2 + in a 2vf}f3 process is allowed by renormalizing 
the GT transition operator with some higher RPA corrections which results in making the 
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matrix elements from Eg. (10. 22 j) non- vanishing. 
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